Noncyclic geomètric changes of quantum states 
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Non-Abelian quantum holonomies, i.e., unitary state changes solely induced by geomètric proper- 
ties of a quantum system, have been much under focus in the physics community as generalizations 
of the Abelian Berry phase. Apart from being a general phenomenon displayed in various subfields 
of quantum physics, the use of holonomies has lately been suggested as a robust technique to obtain 
quantum gates; the building blocks of quantum computers. Non-Abelian holonomies are usually as- 
sociated with cyclic changes of quantum systems, but here we consider a generalization to noncyclic 
cvolutions. We argue that this open-path holonomy can be used to construct quantum gates. We 
also show that a structure of partially defined holonomies emerges from the open-path holonomy. 
This structure has no counterpart in the Abelian setting. We illustrate the general ideas using an 
example that may be accessible to tests in various physical systems. 

PACS numbers: 03.65.Vf, 03.67.Lx 



I. INTRODUCTION 

Berry's discovery [lj of geomètric phase factors associ- 
ated with slowly changing external parameters initiated 
intense studies of geomètric phase effects in quantum me- 
chanics. Wilczek and Zee extended Berry's work by 
pointing out the existence of geomètric effects as a genèric 
feature of quantum adiabatic evolution. They showed 
that the Abelian geomètric phase factors generalize to 
unitary state changes, often referred to as non-Abelian 
quantum holonomies, in the case of Hamiltonians with 
degenerades. Such effects have been shown to occur in 
a variety of systems, ranging from molècules 0, 3 and 
extended systems to atòmic nuclei 0, and fields @ . 
Lately, the interest in holonomies has been refueled due 
to the insight that they can be used to implement fault 
tolerant quantum gates 0, . This has led to work on 
holonomy effects for imp lementations of quantum com- 
putation lli], Il2t lla . Il4 Il5| and quantum information 

EEIïliïMar 

In the aforementioned work, holonomy is associated 
only to loops, i.e., to closed paths of slowly changing 
parameters. But what happens if the path fails to be 
closed? In this work, we address this question and pro- 
pose an approach to noncyclic non-Abelian holonomies. 

Above we have used the language of parameter depen- 
dent Hamiltonians in order to describe the emergence of 
holonomies, where the motion of the eigenspaces of the 
Hamiltonian gives rise to the holonomy. However, the 
Hamiltonian is not a necessary component. We may in- 
stead consider just a subspace moving in the total Hilbcrt 
space of the system. As above, this moving subspace can 
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be realized as an eigenspace of a Hamiltonian, but could 
alternatively be achieved through a sequence of projec- 
tive measurements of observables with the subspace as 
an eigenspace. We use this more general "subspace ap- 
proach" in our definition of open-path holonomy, and ex- 
tend Ref. to the non-Abelian case. However, in order 
to connect with more familiar settings, we also present 
the open-path holonomy in terms of paral·lel transport, 
as well as resulting from adiabatic evolution of parameter 
dependent Hamiltonians. 

Our generalized holonomy contains previous notions, 
such as that of Ref. in the case of cyclically evolving 
Hamiltonians, and that of Ref. 21] for particular paths 
associated with the dynamical invariants of a Hamilto- 
nian. We further demonstrate a concept of open-path 
holonomic quantum gates that may be of use in the 
context of quantum information processing. Finally we 
demonstrate that for some relative orientations of the ini- 
tial and final subspaces of the open path, the holonomy 
is only partially defined. This we call partial holonomy; 
a phenomenon that has no counterpart in the Abelian 
case. 



II. OPEN PATH HOLONOMY 

Consider a smooth curve C in the Grassmann manifold 
G(N; K) [2^, i.e., the set of íí-dimensional subspaces in 
an iV-dimensional Hilbert space. The holonomy for sub- 
spaces should only depend on the properties of this curve. 
There is a natural bijection between the Grassmann man- 
ifold and the collection of projectors of rank K . Thus, 
corresponding to our curve C in the Grassmann manifold, 
we may define a curve P(s) being a family of projectors 
parameterized by s £ [0,1]- Let us now construct the 
intrinsically geomètric quantity 



T = P(1)P(1 - Ss) . . . P{Ss)P(0), 
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where Ss is the step size in a discretization of the curve. 
We are interested in the operator T in the limit of small 
Ss. In order to find an expression for T in this limit we 
let {|afe(s))}jLi be an orthonormal basis of the subspace 
C(s), for each s, and we assume that this family of bases 
is chosen in a smooth way. Note that 

P(s + Ss)P(s)=J2[B(s)}ki\a k (s + Ss)){ai(s)\, (2) 

kl 

where [B(s)]m = S k i + Ss(àk(s)\ai(s)). This allows us to 
rewrite T as 

T = £[B(1 - Ss)B(l - 2Ss) . . . fl(0)]ki|o fc (l)>(oi(0)|. 

kl 

(3) 

We observe that to the first order in Ss, B(s) = 1 + 
SsA(s) = exp (SsA(s)), where 1 is the K x A unit matrix 
and 

[A(s)) kl = (à k (s)\ ai (s)). (4) 
Hence, in the limit Ss — > 0, we obtain 

r = ^[Pe/o^)^] fei |a fc (l))(a / (0)|, (5) 

kl 

where P denotes path ordering. 

A gauge transformation is a change of frames 

\a k (s))^\a' k (s))=J2lU(s)}ik\ai(s)), (6) 
i 

U(s) being a unitary matrix. The set of A-frames, 
i.e., ordered orthonormal A-tuples in an 7V-dimensional 
Hilbert space, forms the Stiefel manifold 

[23- The Stiefel 

manifold can be regarded as a fiber bundle with the 
Grassmannian as base manifold and with the set of A x A 
unitary matrices as fibers. The gauge transformation 
given by Eq. JïïJ) can be seen as a motion along the fiber 
over a point C{s) in the Grassmannian. 

The quantity T is manifestly gauge invariant. On the 
other hand, the matrix Pe^« M s ) ds transforms as 

Peti A{s)ds -» C/t(l)P e i'o A( > s ) ds U(0). (7) 

Hence, the eigenvalues of Pe^» A ( s ) ds are no t gauge in- 
variant 0i|2^| since we may have U(l) =/= U(0). In order 
to deal with this we must somehow find a way to relate 
the initial and final frames. This can be achieved byin- 
troducing the concept of parallel frames |2J, |2jj, |2y, |27j . 

Given a fixed A-frame A — {\a k )} k= i in the subspace 
C a we wish to find a A"-frame B = {\bk)} k= i in the sub- 
space Cb that in some sense is as parallel as possible to 
A. A reasonable approach would be to minimize the fol- 
lowing function over all possible choices of B 

K 

d{ab) = ]rii| afc }-^}n 2 

k=l 

= 2K - 2ReTrM(„4, B), (8) 



where 

[M(A,B)] kl = (a k \bi). (9) 

Thus, in order to minimize D(A, B) we have to maximize 
Re Tr M(A, B) , where it is assumed that B spans over all 
possible A-frames of Cb- We refer to the matrix M(A, B) 
as the overlap matrix. 

Let B — {l^fe}}^! be some arbitrary but fixed K- 
frame of Cb- Every other A-frame B of Cb we may 
write as a unitary transformation of the elements of 
B. All possible overlap matrices can thus be written as 
M(A,B) = M(A,B)V, where V spans over the set of 
unitary A x K matrices. Let M(A, B) = RUm, with R 
positive semi-definite (R > 0) and U m unitary, be a po- 
lar decomposition [23]. If R is positive definite (R > 0), 
then its inverse exists and U m is unique and can be 
constructed as Um = R~ 1 M(A 1 B). 

Note that the positive definiteness of R is a property of 
the pair of subspaces C a and Cb, and not a property of the 
specific choice of frames. In the following we say that two 
subspaces C a and Cb are overlapping if, for any choice of 
frames, the positive part R of the overlap matrix is pos- 
itive definite. One may note that this equivalently could 
be stated as the overlap matrix having A nonzero singu- 
lar vàlues 2§] . In the case when the number of nonzero 
eigenvalues of R is greater than zero but less then A, 
we say the two subspaces are partially overlapping. If 
all the eigenvalues of R are zero, the two subspaces are 
orthogonal. 

In the case when the two subspaces are overlapping 
one can show that the maximum of ReTr(Af(.4, B)V) 
is obtained if we choose V — U' M . Thus, the optimal 
choice of A-frame B is uniquely determined as 

\ik) = Y} u *W)- (10) 

It follows that 

miD{A,B) = 2A - 2Tr^ ' M{A,B)M\A,B) (11) 

with M(A,B) = R. 

An alternative route to find the parallel frame in 
Eq. l(ÏÜ|) is to note that the overlap matrix M{A, B) > 
if and only if B is the parallel frame B. 

If we assume that the initial subspace C(0) and final 
subspace C(l) are overlapping, we can rewrite Eq. (j3J) 
using a final frame {\àk(^))}k=i ^ na t is parallel to the 
initial frame {|afc(0))}^L 1 . This results in 

T = J2[U g ]ki\àk(l))(ai(0)\, (12) 

kl 

where 

U g = U M Peï° Ais)ds . (13) 

Here, Um is the unitary part of the polar decomposition 
of the overlap matrix of the initial frame and the original 
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final frame. Under a gauge transformation of the form 
given by Eq. ©, one can show that the overlap matrix 
transforms as 

M — > M' = U^(0)MU(1). (14) 

This entails that the unitary part of M must transform 
as U M -> U' M = U\0)U M U(1). This fact and Eq. 
cntail that the matrix U g transforms as 

U g ^U' g = UH0)U g U(0). (15) 

Hence, the eigenvalues of U g are gauge invariant and we 
define U g to be the holonomy for subspaces. 

Let us consider some special cases of this holonomy. 
UA{0) = {M0))}f =1 and A{1) = {|a fe (l))}f =1 are two 
parallel frames such that |afc(l)} = \a,k(0)), for all fc, then 
we obtain M(A(0), À(l)) = 1. Hence, U M = 1 and 
U g = Pe-fo A ( s ) ds . Xhis corresponds to the Wilczek-Zee 
holonomy in the case of adiabatic evolution. Further- 
more, when the subspaces are one-dimensional the ma- 
trices reduce to complex numbers. In this case we may 
use Eq. (|13|) to obtain 

U g = e ''arg(<a(0)|a(l)))+/ 1 (à(s)|a(s))ds^ ^g-j 

which fully agrees with the geomètric phase factor in 
Ref. [H. 

Next, we view the holonomy in terms of parallel trans- 
port along the curve C. Intuitively, parallel transport is 
bascd on the notion of transporting a subspace without 
locally rotating it. Assume that we have a family of K- 
frames A(s) = {|ofe(s))}^ 1 parameterized by s S [0, 1]. 
Parallel transport is achieved if and only if A(s + Ss) is 
parallel to A(s), Vs e [0,1). As mentioned above, two 
frames are parallel if and only if their overlap matrix, as 
defined by Eq. JÜJ, is positive definite. The overlap ma- 
trix of the frames A(s) and A(s + Ss) can, to first order 
in Ss, be expressed as 

[M(A(s),A(s + Ss))} M =S kl -Ss[A(s)] kl , (17) 

with A(s) as in Eq. Since A(s) is anti-Hermitian, the 
overlap matrix is positive definite only if A(s) — for all 
s G [0, 1]. Hence, under parallel transport the holonomy 
takes the form U g = Um, where U m is the unitary part 
of the polar decomposition of the overlap matrix between 
the initial frame and the parallel transported final frame. 

Let us now consider adiabatic evolution. Assume H(s) 
is a one-parameter family of Hamiltonians all having 
a degenerate energy eigenspace of dimension K corre- 
sponding to the energy E(s). Furthermore, assume that 
{|a/ c (s))}^L 1 is a basis for the eigenspace. Consider an 
adiabatic change from s = to s = 1 during an clapsc 
of time T. The evolution imposed on a state, initially in 
the degenerate subspace, can be written as 

£7(1,0)P(0) = e - iT ^E{s)ds 

x£[Pe/o A ^ ds } kl \a k (l))( ai (0)l (18) 

kl 



where Í/(1,0) is the evolution operator taking the sys- 
tem from s = to s = 1 and P(0) is the projector 
onto the initial eigenspace. If we assume that the final 
eigenspace is overlapping with the initial eigenspace, we 
may as before consider a final frame that is parallel to 
the initial frame. Using this we may rewrite Eq. i|18[l as 
[7(1, 0)P(0) = e- lT ti E{s)ds T ^ with r ag in Eq The 
first factor of the right-hand side of this equation we rec- 
ognize as the dynamical phase factor, while the second 
contains the open-path holonomy. 

The total action of T in Eq. I|12fl can be decom- 
posed into two parts. One part is given by the par- 
tial isometry T — X^fc=i l^fc(l))(«fc(0)|, which maps the 
initial frame to its parallel frame. The second part is 
R = J2k.i=i{Ug}ki\àk(í)){ài{l)\, which is a unitary trans- 
formation on the final subspace given by the holonomy. 
This decomposition of T provides an understanding of 
how the holonomy should behave under a gauge trans- 
formation. In order for the unitary transformation on the 
final subspace to be independent of gauge, the holonomy 
must transform unitarily, as displayed in Eq. i|15|) . 

In the language of quantum computation, one may say 
that we choose to let the parallel frame {|afc(l))}^L 1 in 
the final space correspond to the computational basis 
{la/^O)}}^! in the initial space. The holonomy then 
appears as the resulting operation with respect to these 
choices of computational bases. An aspect of this is 
that the computational basis becomes path-dependent. 
One might, as an example, consider a sequence of open- 
path holonomic implementations of operations. If this 
sequence happens to join into a cyclic path, it might be 
the case that the initial computational basis does not co- 
incide with the final computational basis, although they 
span the same subspace. 



III. PHYSICAL EXAMPLE 

In order to illustrate the concept of open-path holon- 
omy, as well as to provide an explícit example of an open- 
path holonomic implementation of a single qubit gate, we 
now consider a specific model system. This model was 
first examined in connection to non-Abelian holonomies 
in Ref. [2^| and would be accessible to tests in vari- 
ous physical systems, such as ion traps [ÏJ, |2íj, àtoms 
[Ï3| , superconducting nanocircuits [Ï4]| , and semiconduc- 
tor nanostructures [ïj| . The Hamiltonian of the system 
reads 

H(s) =wb(a)|e)(0|+«i(s)|e>(l|+a;,(s)|e>(a|+H.c., (19) 

where |0),|l),|a), and |e) are orthonormal, and íüq(s), 
üJi(s), and u> a (s) are tunablc coupling parameters. We 
assume that the parameters combine to a real vector 
(Lüo(s),u>i(s),uJa(s)) of unit lcngth. Thus the parameter 
space forms a unit 2-sphere, which we may parametrize 
using the polar angle 9 and the azimuthal angle ip. The 
Hamiltonian H(s) has a doubly degenerate zero-energy 
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eigenspace, which is spanned by the eigenstates 



IV. PARTIAL HOLONOMY 



\Di(s)) — cos 9(s) cos ip(s)\0) + cos 8(s)sinip(s)\l) 
— sm.9(s)\a), 

\D 2 (s)) = -sin^(s)|0)+cos^>)|l}, (20) 

where 6(s) G [0, n] and ip(s) e [0, 2tt). In this context, 
the states \Di(s)) and \D2(s)) are often referred to as 
"dark states" . 

Let us now assume that the parameter s is changed 
slowly enough for the evolution to be adiabatic. Further, 
let (9(0)^(0)) = (0,0) and (0(1), </j(1)) = (0i,<^i). The 
overlap matrix between the initial frame .4(0) = {|0), |1)} 
and the final frame ,4(1) = {|-D X (1)), |L> 2 (1))} is 



M(A(0),A(l)) = 



cos t'í cos (pi — sint^i 



l cosf^siní^! cos ipi J ' 
The unitary part of the overlap matrix is 

- siní^jA 



U M 



COS Oi COS pi 

| cos 9± | 

COS 8-\ Sill 



(21) 



(22) 



which holds under the assumption that 9\ ^ -n/2. If 
we furthermore assume that < 9\ < ir/2, we may write 
U m = e~ lVlCr y , where <j v is the y-component of the Stan- 
dard Pauli matrices. For the frame in Eq. Q20|l . we obtain 
A(s) — i cos 9 (s) (p (s)(T y , yielding 



U — e - itT v{<Pi-fo cos 0{s)^(s)ds) 

— e - ier v{<Pi-fc cosed v) _ p-ï<r»7 



(23) 



where the quantity 7 equals the sòlid angle swept by the 
geodèsic closure of the curve C on the parameter sphere. 
If we instead assume that 7r/2 < 9\ < n, we obtain 



U 



M 



cosi^i -sin^i 
-sin(/?i cost/?i 



(24) 



which can be written as Um = e~ iVl ' Ty (—a z ), where a z 
is the z-componcnt of the Standard Pauli matrices. In 
this case the holonomy takes the form 



-i<Pl<ry \ e ~i<r V Se cosedí P^ 



(25) 



Due to the fact that the different components of the Pauli 
matrices do not commute, the holonomy is no longer de- 
termined by the sòlid angle swept by the geodèsic clo- 
sure of the curve C on the parameter sphere. In the 
first case the holonomy had an Abelian structure (in the 
sense of Ref. [13) due to the fact that [A(s),A(s')] = 
for any s,s' 6 [0,1] and [U M , eti A ( s ~> ds ] = 0. The lat- 
ter is not fulfilled in the second case, were the holonomy 
is truly non-Abelian. Hence, for this physical example 
open paths seems to be a necessary prerequisite in order 
to achieve truly non-Abelian holonomies. 



So far we have assumed that the initial and final sub- 
spaces of the open path are overlapping. In the spe- 
cial case of a one-dimensional subspace there are two 
cases, either the subspaces are overlapping, or they are 
orthogonal. As a conscqucncc the holonomy either ex- 
ists uniquely, or is undefined. In the non-Abelian case, 
however, there is an additional case, namely that the sub- 
spaces are partially overlapping. In this case the holon- 
omy is only partially determined. When the two sub- 
spaces are partially overlapping the positive part R of the 
overlap matrix is not invertible, no matter the choice of 
frames. However, we may use the Moore-Penrose pseudo 
inverse (MP-inverse) |2g. Since R is a positive semi- 
definite matrix, its MP-inverse _R e can be calculated by 
inverting the nonzero eigenvalues in its spectral decom- 
position. The matrix U m can now be defined as the 
partial isometry U m = R e M. This results in a partial 
isometry R e M Pe^o A ^ ds that we shall call a partial 
holonomy. 

Let us examine how the partial holonomy behaves un- 
der a gauge transformation. The overlap matrix between 
the initial and final subspaces transforms as in Eq. H14[) . 
It follows that 



R' =VM'M 



'71 í't 



C/ t (0)Mí7(l)í7 t (l)M t í7(0) ( 26 ) 
C/t (0)JH7(0), 



where R = We note the following property 

of the MP-inverse. Suppose that U and V are unitary 
matrices. Then, for any matrix X, we have (see p. 434 
in Ref. [H) 

(UXVf = v^x e u^. 

Thus, 

U' M =(R') e M' 



(27) 



17* (0)JH7(0) ] U\0)MU{\) (28) 



=[/ f (0)fí e M[/(l), 

which is precisely the way R~ X M transforms if R is 
invertible. Hence, the transformation of U m takes the 
same form independently of whether or not the matrix 
R is invertible. Moreover, the path ordered part of the 
holonomy, Pe^o M s ) ds ; always constitutes a unitary ma- 
trix that transforms according to Eq. Q . Thus the par- 
tial holonomy transforms unitarily just as the holonomy, 
as displayed in Eq. fTÏÏf. 

As an example of a partial holonomy, let us revisit the 
previous model system, now assuming that 9± = tt/2. 
The overlap matrix M(A(0),A(1)) reduces to 



M(A(0),A(1)) 



Q, Q = 



o A 
1, 



(29) 
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and R = e ~ lVirTy Qe l·ipirTy , which happens to be a one- 
dimcnsional projector, and thus equal to its own MP- 
inverse. Consequently, Um — M(A(0),A(1)) and the 
partial holonomy becomes 

U = g-Wifs Q e i lo cos 8{s)ip(s)ds<r y ^ ^gg-j 

One may note that the existence of a loop (in this case 
the equator 9i = tt/2 of the parameter sphere) along 
which the holonomy is not fully defined is a topologi- 
cally enforced prerequisite for the discontinuous transi- 
tion between the Abelian and non-Abelian character of 
the holonomy, that wc have found in this examplc. 

V. CONCLUSIONS 

To summarize, we consider subspaces moving in the 
Hilbert space of a quantum system, and the concomitant 
unitary transformation associated with the geometry of 
the traversed path; the non-Abelian quantum holonomy. 
The Standard non-Abelian quantum holonomy is defined 
for closed paths of such subspace motions, while we con- 
sider an open-path generalization. Due to the opcnness 



of the path, the initial and final subspaces do not coincide 
in general. In order to "extract" the unitary transforma- 
tion on the final subspace, i.e., the holonomy, we use a 
concept of parallelity in order to decide which basis in 
the final subspace corresponds to the basis in the initial 
subspace. Under suitable conditions on the relative orien- 
tation between the initial and final subspaces, this proce- 
dure results in a uniquely defined non-Abelian quantum 
holonomy for open paths. This enables the construction 
of quantum gates in the open-path setting, where the ac- 
tion of the gates is given by the proposed holonomy. The 
idea of open-path holonomic gates may be useful when 
analyzing noncyclic errors 30, 31] of Standard implcmcn- 
tations of holonomic quantum computation. We finally 
point out the existence of partially defined holonomies, 
which has no counterpart in the Abelian case. 
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tum Computation at DAMTP, Cambridge, for hospital- 
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Swedish Research Council. 
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